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SUMMARY 


Using superelements (substructures), modal and frequency response 
analysis was performed for a large model of the Advanced Attack Helicopter 
(AAH) developed for the U. S. Army. Whiffletree concept was employed so 
that the residua] structure along with the various super elements could be 
represented as beam-like structures for economical and accurate dynamic 
analysis. A very large DMAP alter to the rigid format was developed so that 
the modal analysis, the frequency response, and the strain energy in each 
component could be computed in the same run. 


INTRODUCTION 


The helicopter model shown in figure 1 consists of 13 substructures 
having 1000 grid points and 4000 structural elements. The dynamic analysis 
was performed using the residual structure, which represented a beam-like 
structure along the centerline of the aircraft at its nominal elastic axis, with 
appendages in the form of other beams. The residual structure was connected 
to the super elements by rigid elements in a sort of whiffle tree'* form. This 
is an interesting, newly developed concept in that it allows dynamic analysis 
with the efficiency of a "stick" model, resulting in significant cost saving 
while at the same time retaining more involved effects accounted for in the 
detailed finite element model. 
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In order that normal mode analysis as well as frequency response 
computations could be executed in the same run, a special alter to the rigid 
format was developed. The alter was extended to further include calculations 
of strain energies in each component, so that the critical parts of the structure 
for each model could be easily identified. To determine the most sensitive 
mode to the rotor excitation, rotor impedance factors were calculated using 
eigenvectors generated from NASTRAN run. 

The structural contributions due to the different subsystems, designed 
by other companies, were incorporated into the main helicopter model by 
employing the NASTRAN general elements (GENEL), To better evaluate the 
interaction between the basic structure and the subsystem, transient response 
analysis due to gun firing forces and frequency response analysis due to the 
rotor excitation were conducted by utilizing the results computed by the 
normal mode analysis. 

Separate set of values obtained from the analytical approach were found 
to be in very good agreement. Harmonic analysis is being currently per- 
formed using the time history output from die transient response analysis 
to determine tlie frequency content. 


DESCRIPTION OF THE MODEL 


The Advanced Attack Helicopter consists of the fuselage, the wing and 
the vertical tail, all composed of typical skin-stringer construction. The 
horizontal stabilizer is a semimonocoque structure utilizing beaded skin 
cons truction. The complete finite element model is shown in figure 1* This 
consists primarily of simpler elements in NASTRAN library. This model 
was divided into 10 primary super elements and three secondary (image) 
super elements. Advantages of substructuring and its limitations were fully 
considered in selecting this particular scheme of dividing up the model. 
Figures Z and 3 show a more detailed definition of a typical superelement 
representation. Table I provides the element breakdown for each of the 
super elements in the model. As can be seen, this is definitely not a small 
model for dynamics. For the convenience of the users, table II shows the 
various set sizes for each super element. 

The dynamic analysis is performed using the a-set points, which in the 
superelement concept is defined as the residual structure or superelement 0, 
This structural model is a beam-like structure along the centerline of the 
vehicle with appendages in the form of other beams and/or GENEL elements. 
This residual structure is not a normal structure consisting of a-set points 
but is connected to the other superelements by means of a unique set of rigid 
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elements. This concept is given the name of ^‘whiffle tree, “ further described 
in the next section. 


Wr-IIFFLETREE CONCEPT 


Dynamics analysis involves some form of technique to reduce the size of 
problem from the sl.^tics analysis model. Most people are very familiar with 
the “stick"or beam or El/ GJ models and even today are quite content using 
those models. With NASTRAN users, Guyan reduction has become quite a 
popular technique for achieving the smaller model. Whiffletree concept uses 
rigid elements or, more appropriately, MFCs for the same goal; namely, 
tlie smaller model. MFC or the multipoint constraints can be used such that 
the behavior of several grid/mass points is represented by one single grid/ 
mass or a scalar point. 

Figure 4 shows a typical whiffletree arrangement for one of the bulkheads 
in tlie AAH model. As can be seen, tlie point in tlie middle (does not have to 
be located there) is ’’connected^* to the more important grid points on the bulk- 
head using a general rigid element RBE3. This single poim: now represents 
tlie average behavior of tliat complete bulkhead without any matrix reduction. 
Continuing in this fashion the end xDroduct is a pseudo stick model tliat would 
predict the dynamic characteristics of the complete structure. Figure 5 
shows this model. 

These points were all left in the residual or boundary structure and thus 
tliey were in the a-set of each superelement. As a matter of fact, the residual 
structure consists of this kind of points only. Because of tliis technique it was 
possible to keep the sizie of a-set points well witlain control for high efficiency 
while retaining in the solution the more complex behavior of the total structure. 
Certainly MCEl and MCEZ modules are not that inexpensive but, compared 
to tlie exponential cost increase of SMPl and SMFZ modules, this approach 
was found to be more cost-effective. 


ALTERS FOR CONVENIENCE 


One of the major tasks undertaken to complete this project consisted of 
formulating an Alter to the Rigid Format so that the following tliree analyses 
could be performed in the same z’un: 


(a) Normal modes analysis witli frequency 
response, 
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(b) Rotor impedance factors computation, 

(c) Modal strain energy distribution. 

All three alters are merged into one large alter to make it a production 
tool. ’Altliough it is not easy to keep up tlie alter compatible with the newer 
versions of MSC/NASTRAN, successful transition was made from CDC 6600 
Version 32 to IBM Version 38; and, presently, efforts are in progress for 
adopting the alter to the latest IBM Version 46. 


Natural Frequencies and Modes 

As shown in table II, the total dynamic degrees of freedom (a-set) add 
up to tlie matrix size of 177. Grid points witli dynamic degrees of freedom 
(a-set) are shown in table II. The modes are identified in somewhat arbitrary 
manner by looking at the deformed shape of tlie structure in a particular mode 
in conjunction with tlie computed mode shapes and strain energy distribution. 
The isometric views of the deformed structure are shown in figures 6 and 7, 


Rotor Impedance Factors 

The rotor impedances are calculated from tlie relation 


1 . 2 
7 : 75 a 

rotor impedcuice i 


where is the eigenvector at the main rotor head in tlie i^^ direction, and 
are plotted in figure 8 as versus the significant modes in their relative 
order of modal strengtli. The modes with tlie lowest impedance (i. e, , tlie 
highest inverse ratio) are most sensitive to rotor excitation. 


Modal Strain Energy Distribution 

The MSC/NASTRAN computer program witli tliis alter lias tlie capability 
to compute and print out strain energy distribution in each flexible mode. This 
is very helpful in identifying tlie critical parts of tlie structure for each mode. 

The program computes tlie strain energy in percentage form in each 
superelement for the desired modes. 
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The results are shown by MATRIX SPT in figure 9- The column number 
identifies the mode number (e. g, , col 7 means mode 7) and the rows 1 through 
10 denote the corresponding superelements 1 through 10* The strain energy 
in Eiipereleinent number zero is computed separately under MATRIX RPT; 
however, it is included in figure 9 for convenience. 


FREQUENCY AND TRANSIENT RESPONSE 


Frequency Response 

Alternating aerodynamic forces acting on rotor blades and on the fuselage 
and nonrotating parts of the vehicle are tlie major source of vibrations* Varia- 
tions in these forces are periodic and all tlie steady alternating force inputs 
to the rotor hub occur in even multiples of the rotating speed such as i/rev 
(IP), 2/rev (2P), 3/rev (3P), etc. However, only alternating forces and 
moments which are integral multiples of the number of blades are transmitted 
to the rotor hub. This helicopter has four-bladed rotor, therefore 4/ rev, 

8/ rev, etc* , are the only input to the hub* The 4/rev (4P) excitation forces 
which are the major contributor to the input at the hub were considered for 
tliis analysis. These excitation forces are obtained from DART (Dynamic 
Analysis Research Tool) program for various forward speeds of the helicopter 
and are multiplied by 4/rev response loads per unit excitation at rotor hub 
{computed by NASTRAN PROGRAM) to obtain total response at crew, stabilized 
sight and various other desired locations. 

The effect of rotor speed variation on pilot and copilot station vibrations 
assuming no change in rotor forces is shown in figures 10 and 11. 


Transient Response 

Transient response due to gun firing is computed at several locations, 
such as crew stations, stabilized sight and gun center of gravity* The input 
impulse time history is shown in figure 12, and a typical response time 
history in figure 13. 

CONCLUSIONS 


A large model for dynamic analysis has been successfully used for the 
AAH project. Once again it is clear tliat NASTRAN is capable of solving a 
rather complex analysis scheme in a production manner provided proper 
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resources are put into planning and writing some interesting DMAPs, The 
concept of MFCs and/or rigid elements is a very powerful tool that seemfj to 
provide new answers to several old problems. 
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TABLE 1 - VARIOUS ELEMENTS IN THE MODEL 
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TABLE II - VARIOUS MATRICES IN THE MODEL 


SLiperclemoiit 

No, 

Matrix 

Siiie 

1 

— 

2 

3 

4 

5 

!> 

7 

8 

1 

10 

0 

G Y i d 
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140 

37 


50 

7 

27 

27 
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1266 

840 

222 

222 

300 

42 

162 

162 

870 

870 

ni - 3 ol 
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42 

54 

30 

30 

60 

0 

68 

68 

134 

36 

f-sol 

1210 

774 

511 

95 

95 

240 

36 

68 ^ 

68 

430 

648 

o- s ct 
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^97 

59 j 

59 

136 , 

24 1 

56 

56 

382 

471 
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i 
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36 

36 

54 

12 
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1 

12 
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1 
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Figure 1. Advanced Attack Helicopter Model 
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Figure 2. Advanced Attack Helicopter Tail Cone (Super element 3) Model 
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Figure 3. Advanced Attack Helicopter Model 
Empennage (Superelement 10) 
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Figure 4. Typical Whiffletree Arrangement 
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Figvire f>. First: Vorlicti-l Mode 
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Figure 8. Rotor Impedance Factors 
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ADVANCED ATTACK HELICOPTER, 
MODEL 41, 4 ROCKET PODS 


APRIL 2S, 1978 


SAME AS MOD. 40 EXCEPT WITH 4 ROCKET STORES 

MATRIXSPT (GTNO name 101) IS A REAL 50 COLUMN X 1 1 ROW RECTANG MATRI X. 
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Figure 9. Typical Strain Energy Output 




Figure 10. Frequency Response at Pilot Seat 
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Figure 11. Frequency Response at Copilot Seat 
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